
TITLE 


AUTHOR 


COMPANY 


DATE 



SOURCE LANGUAGE 



PROGRAM LIBRARY 


8-568 

CFI - CONTINUED FRACTION INVERSION 

Andres T. Siy 

Capitol Institute of Technology 
Kensington, Maryland 

July 14, 1972 

8K FORTRAN 


Although this program has been tested by the contributor, no warranty, express or implied, is made by the contributor. 
Digital Equipment Computer Users Society or Digital Equipment Corporation as to the accuracy or functioning of the 
program or related program material, and no responsibility is assumed by these parties in connection therewith. 






DECUS ftogram Library Write-up 


DECUS NO. 8-568 


Sly 


SUB 8K-# 


Jan. 10, 1972 


( 1 ) NAME: CPI Continued Fraction Inversion 

—tasesssbo convert a continued fraction into a 
rational function. 

(2) CALLING SEQUENCE : 

CALL CFI (N,H,B,A) 
where : 

Input Data : 

N a one-half the no. of given H»e (SEE ( 4 ) & (7) ) 
H(I),I - 1 , 2 , ..., 2 *N (SEE ( 4 ) & (7) ) 

Return Data : 

B(I),I * 1,2, 00 ., H coefficients in numerator polyno. 
A(I),I a 1,2.(N+1) " ” denominator " , 

( 3 ) ERROR RETURN : 

-none- 

(4) SPECIAL CONSIDERATION : 

See (7) 

(5) SUBPROGRAM CALLED : 

( 6 ) STORAGE REQUIRED : 

pages (octal) 

( 7 ) ALGORITHM & REFERENCES : 

Given: G(s) » v *___:--——— 
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Algorithm : 


1 for i a 1 j> k 31 , »»* g 1 


Step 1: define a ilc c q £ 0r £ greater than*k?st ,2,.. ,n*1 


to ik “ 0 if i is greater than k 

Step 2: For each k «2,3, ...» n*i compute for i *2,3,«.,k 

b ik 3 h 2(n~k)+4 a ±-1,k-1 + b i,k-t 

^ik 3 h £(n~k)*3 b lk * *1^-1 

» 

Step 3? compute, 

h i ’ b l*1 ,u 

®i * ®i»1 ,u 
for i * 1,2, oov>, n 


References : 

T. C,To Chen, - ; >"A Formula and an Algorithm for 

Continued Fraction Inversion", Proc. XEEE(letters) 
vol. No. pp. 1780-1781, Oct. 1969. 

2. C.S?. Chen, L»S. Shieh, "Continued Fraction Inversion 

by South*s Algorithm ", IEEE Trans Circuit Theory 
vol. CT-lSj, pp. 197-202, May 1969. 


(8) LISTING : 

See attach* 

(9) SAMPLE : 

See attach. 
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**** " CONTINUED FRACTION INVERSION " 

**** RESULTS IS A RATIONAL FUNCTION 

**** CBC1>*S**CN-1> + ••• + 8CN)>/CS**N + AC 1>*S**CN-l) + 
••••+ACN)) 

JAN* 10*72 
N IS ATMOST 10 
SUBROUTINE CFI CN#H/9*A> 

DIMENSION HC20)#B C10 ) #ACl 1)>A1C11>11)#B1C11#11 > 

N1»N+1 

DO 10 J«1#N1 
A1C1,J)»1* 

Ji=J+i 

DO 10 I»J1#N1 
A1CI*J)*0* 

B1CI#J)«0. 

10 CONTINUE 

DO 80 K=2#N1 
DO 20 I»2*K 
NK=2*CN-K*2> 

BlCI#K)«HCNK)*AICI-t#K-l> + B1CI>K-1> 

A1 CI jK)»HCNK-1 )*B1CI J |K)+A1CI#K“1> 

20 CONTINUE 

DO 33 1=1#N 
11=1 + 1 

BC1)*8lCI l*Nl) 

AC 11)®A1C1I#N1) 

33 CONTINUE 

AC1>«l. 

RETURN 

END 
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C ***♦ SAMPLE f 1 

C *** HR ■ DUMMY ALPHANUMERIC DATA MUST 3E A 
DIMENSION KC14>*BC7>#ACR) 

3 nEAD<f»*5> nn#N 

5 KOHMATCAl#110) 

IF CN> 100, 103, 10 
13 NX»2*N 

DO 20 I»1,NN 
20 KEADC2,15) HCI) 

15 FORMATCE10*4> 

C 

CALL CF1(X*H*B,A) 

*>KITEC2,35)KR, 

DO 33 I *>1 *X 

30 WniTE<2,15> BCD# 

FINI 

35 FORMAT C//9X, AD 

N!-N+1 

DO 40 I«l#Nl 

43 it/RITEC2#l 5) AC D, 

FIX I 
GO TO 3 
100 STOP 

END 
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*20002+01 


•1000E+01 

•2000E+31 



• 6900E+01 

•2400E+02 

•1000E+31 

•1830E+02 

•2400E+02 


•1200E+02 

•2400E+03 •7200E+03 

•1003E+01 

•7200E+02 

•6303E+03 .7200E+03 


•2000E+02 

• 1200E+04 *1512E+05 

• 10PJE+01 

•2030E+03 

•5400E+34 *3528E+05 


•3033E+02 

•4200E+04 *141 l'E+06 

•1000E+31 

•4530E+03 

•2940E+05 *5645E+06 


•6030E+01 

•3000E+00 

•1000E+01 

•6030E+01 

•3030E+00 


•4000E+01 

•0000E+00 

•1000E+01 

•4000E+01 

•0000E+00 


•1500E+01 

*4000 E+0 lCVCC)— \ + 

• 1300E+01 

•8003E+31 

•1200E+02 V * 


•3000E+00 

•0000E+00 

.1000E+01 

•0000E+00 

•0000E+00 


(V* w^iry-) 
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